Probability density function characterization of multipartite entanglement 



P. Facchi,^'^ G. Florio,^'^ and S. Pascazio^'^ 

' Dipartimento di Matematica, Universita di Bari, 1-70125 Bari, Italy 
^INFN, Sezione di Bari, 1-70126 Bari, Italy 
Dipartimento di Fisica, Universita di Bari, 1-70126 Bari, Italy 
(Dated: February 1, 2008) 

We propose a method to characterize and quantify muhipartite entanglement for pure states. The 
method hinges upon the study of the probabihty density function of bipartite entanglement and is 
tested on an ensemble of qubits in a variety of situations. This characterization is also compared to 
several measures of multipartite entanglement. 
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I. INTRODUCTION 

Entanglement is one of the most intriguing features 
of quantum mechanics. Although it is widely used in 
quantum communication and information processing and 
plays a key role in quantum computation, it is not fully 
understood. It is deeply rooted into the linearity of quan- 
tum theory and in the superposition principle and basi- 
cally consists (for pure states) in the impossibility of fac- 
torizing the state of the total system in terms of states 
of its constituents. 

The quantification of entanglement is an open and 
challenging problem. It is possible to give a good def- 
inition of bipartite entanglement in terms of the von 
Neumann entropy and the entanglement of formation 
[H . The problem of defining multipartite entanglement is 
more difficult and no unique definition exists: differ- 
ent measures capture in general different aspects of the 
problem Attempts to quantify the degree of quan- 
tum entanglement are usually formulated in terms of its 
behavior under local operations/actions that can be per- 
formed on different (possibly remote) parts of the total 
system. Some recent work has focused on clarifying the 
dependence of entanglement on disorder and its interplay 
with chaos [j, [al , or its behavior across a phase transition 

m- 

The work described here is motivated by the observa- 
tion that as the size of the system increases, the num- 
ber of measures (i.e. real numbers) needed to quantify 
multipartite entanglement grows exponentially. A good 
definition of multipartite entanglement should therefore 
hinge upon some statistical information about the sys- 
tem. We shall look at the distribution of the purity of a 
subsystem over all possible bipartitions of the total sys- 
tem. As a characterization of multipartite entanglement 
we will not take a single real number, but rather a whole 
function: the probability density of bipartite entangle- 
ment between two parts of the total system. The idea 
that complicated phenomena cannot be "summarized" 
in a single (or a few) number(s) stems from studies on 
complex systems Q and has been considered also in the 
context of quantum entanglement In a few words, 
we expect that multipartite entanglement be large when 



bipartite entanglement is large and does not depend on 
the bipartition, namely when its probability density is a 
narrow function centered at a large value. This character- 
ization of entanglement will be tested on several classes 
of states and will be compared with several measures of 
multipartite entanglement. 



II. THE SYSTEM 

We shall focus on a collection of n qubits. The di- 
mension of the Hilbert space is iV = 2" and the two 
partitions A and B are made up of ua and ub spins 
(jiA+nB — n), respectively, where the total Hilbert space 
reads H = Ha ® Hb and the Hilbert spaces Ha and Hb 
have dimensions Na = 2"-* and Nb = 2"^, respectively 
{NaNb = N). We shall consider only pure states 



(1) 



where \k) = IJa) 'Si \Ib), with a bijection between k and 
UaJb), < Ja < Na - 1 a.nd < Ib < Nb - 1. As a 
measure of bipartite entanglement between A and B we 
consider the participation number 



Nai 
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7I-AB Tr^ p^, PA = Trs p, 



(2) 



where p = \ip){'ijj\, and Tr^ (Trs) is the partial trace 
over the degrees of freedom of subsystem A {B). Nab 
can be viewed as the relevant number of terms in 
the Schmidt decomposition of IV") The quantity 

riAB = log2 Nab represents the effective number of en- 
tangled spins. Clearly, for a completely separable state, 
Tr^ p\ = 1 for all possible bipartitions, yielding Nab = 1 
and uab = 0. In this sense the participation number 
can distinguish between entangled and separable states. 
Moreover ttab is directly related to the linear entropy 
Sl = 1 — ttabi that is an entanglement monotone, i.e. it 
is non increasing under local operations and classical 
communication. In general, the quantity Nab will de- 
pend on the bipartition, as in general entanglement will 
be distributed in a different way among all possible bi- 
partitions. Therefore, its distribution p(Nab) will yield 
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information about multipartite entanglement: its mean 
will be a measure of the amount of entanglement in the 
system, while its variance will measure how well such en- 
tanglement is distributed, a smaller variance correspond- 
ing to a higher insensitivity to the particular choice of 
the partition. 

We will show that for a large class of pure states, sta- 
tistically sampled over the unit sphere, p{Nab) is very 
narrow and has a very weak dependence on the biparti- 
tion: thus entanglement is uniformly distributed among 
all possible bipartitions. Moreover, p{Nab) will be cen- 
tered at a large value. These are both signatures of a 
very high degree of multipartite entanglement. 

By plugging ([T|) into ^ one gets 



Na-1 Nb-1 

E 



E 



ZjlZj'lZj/l'Zjli 



(3) 



We note that ttab — T^'caPa — Ttb p% and 1/Na < 
Ttap\ < 1, with the minimum (maximum) value at- 
tained for a completely mixed (pure) state pA- There- 
fore, 



1 < Nab = Nba < min{iVA, Nb}- 



(4) 



A larger value of Nab corresponds to a more entangled 
bipartition {A,B), the maximum value being attainable 
for a balanced bipartition, i.e. when ua = [n/2] (and 
n-B — [("^ + l)/2]), where [x] is the integer part of the real 
X, that is the largest integer not exceeding x, and the 
maximum possible entanglement is Nab = Na = 2"-* = 
VN {— N/2) for an even (odd) number of qubits. As 
anticipated, as a characterization of multipartite entan- 
glement we will consider the distribution of Nab over all 
possible balanced bipartitions. 



III. MEASURING MULTIPARTITE 
ENTANGLEMENT: SOME EXAMPLES 

Let us illustrate this approach on the simplest non- 
trivial situation, that of three entangled qubits. If the 
pure state is fully factorized, say 



W = \k) 



(5) 



for a given < /c < 7, then the reduced density matrix 
PA of every qubit is a pure state, whence 



P{Nab) = Snai 



(6) 



there is no entanglement. On the other hand, for a max- 
imally entangled state 
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|^) = ^(10002) + IIII2)), 



(7) 



one gets a completely mixed state for every partition, 
namely = I2/2 and thus 



P{Nab) = Snae 



(8) 



with maximum average and zero variance: there is maxi- 
mum multipartite entanglement, fully distributed among 
the three qubits. The above probability distributions 
should be compared with an intermediate case like 
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IV) = ^(10002) 



IIIO2)), 



(9) 



where the first couple of qubits are maximally entangled 
(Bell state) while the third one is completely factorized. 
In such situation one gets pi = p2 = I2/2, while pa = 
|1)(1|, whence 



P{Nab) = '^jv.iB.l/S + 2^Ar^B,2/3. 



(10) 



This simple application discloses the rationale behind the 
quantity p{Nab) as a measure of multipartite entangle- 
ment. 

When the system becomes larger, the natural exten- 
sion is towards larger (balanced) bipartitions. We stress 
that, besides the comment that follows Eq. the use 
of balanced bipartitions is simply motivated by the fact 
that, in the thcrmodynamical limit, the unbalanced ones 
give a small contribution, from the statistical point of 
view: this can be easily understood if one considers that 
for n large and ua <IC n the binomial coefficients 



n 
n/2 



> 



n 

UA 



(11) 



so that our characterization of multipartite entanglement 
will be largely dominated by balanced bipartitions. No- 
tice also that very unbalanced bipartitions of large sys- 
tems yield negligible average entanglement [l2| [lol ■ For 
all these reasons, if one considers the distribution over 
all bipartitions, the contribution from the balanced bi- 
partitions will dominate due to PT|) . By contrast, if only 
unbalanced bipartitions are considered the results will be 
in general very different. 

It is interesting to study the features of the characteri- 
zation of entanglement proposed in Sec. [Ill when applied 
to particular classes of states. For the GHZ states 
we find 



Nab{GRZ) = 2 



(12) 



for all possible bipartitions (both balanced and unbal- 
anced) and for an arbitrary number of qubits. Clearly, 
the width of the distribution is 0, i.e. p{Nab) = Snab,2- 
For the W states [111 we obtain 



Nab{W) 



(13) 



This value depends only on the relative size of the two 
partitions, i.e. also in this case the width of the distri- 
bution of bipartite entanglement is 0. Notice that, if n 
is even, Nab{^) ~ 2 for balanced bipartitions (and in 
this case a discrimination between W and GHZ states 
would require the analysis of unbalanced bipartitions). 
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Moreover, in the large n limit NabO^) — 2 also for n 
odd. 

These results indicate that, for n large, the amount of 
(multipartite) entanglement is limited both for GHZ and 
W states. These states essentially share the same amount 
of entanglement when n is large. They can be distin- 
guished only by considering less relevant (from the sta- 
tistical point of view) bipartitions. Moreover, for n large, 
-^ab(W) ^ 1 for balanced bipartitions. This means that 
also in the thcrmodynamical limit the W states retain 
some entanglement. 



IV. TYPICAL STATES 

Let us now study the typical form of our characteriza- 
tion of multipartite entanglement p{Nab) for a very large 
class of pure states of the form ^ , sampled according to 
a given statistical law. Several features of these random 
states arc already known in the literature [T^, , but 
we shall focus on those quantities that are relevant for 
our purpose. We write 



N-l 



(14) 



k=0 



where (jjk are independent random variables with expec- 
tation 



E[e"l"'] = 



(15) 



and r = (ri,...,rjv) is a random point with a given 
symmetric distribution p{r) on the hypersphere 5^"^ = 
{r G R^jr^ = 1}. The features of these random states 
are readily evaluated: one first splits ttab in two parts 



ttab = Xab + Mab , 



(16) 



where 



(17) 

= E'E''>''' + EE''-M''+E4' (18) 

j-j' I j l-l' 3,1 

with j,/ = 0, . . . , 7V^— 1, 1,1' ~ Q, . . . , Nb — 1; and primes 
banning equal indices in the sums. 

We note that the expectation value i?[r|;] = 0{1/N), 

thus Xab and Mab are sums of at most N'^ terms of 
order 1/N'^. By the central limit theorem, for large N, 
-KAB tends to a Gaussian random variable with mean and 
variance 



II AB = E[kab\, 

2 irr 2 1 2 



(19) 



respectively, namely it is distributed as 



/(ttab) 



1 / [T^AB-fJ.AB) 1 



From E[Xab] = and the independence between phases 
(pk and moduli we get 

liAB - E[Mab] = N^Na + Nb~ 2)E[rlri] + NE[rt] 

(21) 

and 



<B^E[XlB]+E[MlB]~^^lB, 



where 



E[Xab] 



2N{Na-1){Ne 



l)E[rlrjrlrl] 



(22) 



(23) 



and 



E[Mls] = N{Na + Nb-2) 

X [{Na + NB)iN - 4) - 2(iV - 5)]E[rlrlrlrl] 
+2N{Na + Nb-2){N + 2Na + 2Nb - 8)E[rlrlri] 
+N{N + 2Na + 2Nb - 
+iN{NA + NB-2)E[rir: 



5)E[rfrl] 



NE[rl] 



(24) 



where we used ii'[r"r2 ^grl^ 



,4j - E[rfr^^r]Ti] witht,j,l,k 
all distinct. Notice that the above results do not depend 
on the particular distribution of 0^, as far as the condi- 
tion is satisfied (otherwise the analysis is still valid, 
but Eqs. (PT|) - ([M1) become more involved). Our results 
particularize for the case of a typical pure state ([1]), sam- 
pled according to the unitarily invariant Haar measure, 
where each Zk € C is chosen from an ensemble that is 
uniformly distributed over the projective Hilbert space 
X^fcl-^feP = 1- such a case, in (fH|) . (f)k S [0, 27r] are 
independent uniformly distributed random variables and 
r = (r 1, r^r) is a random point uniformly distributed 
on the hypersphere S^~^, with distribution function 



■)N 



p{r) 



tN/2 



(25) 



the prefactor being twice the inverse area of the hyper- 
octant {ri > 0}, with r(.T) the Gamma function. 

The explicit expressions of pT|) -p4 |) can be computed 
through (|25p , recovering the values of mean and variance 
obtained by different approaches [1, [l^, [l^. However 
one can easily estimate them for large N by the following 
reasoning. For large N the marginal distributions of the 
amplitudes rk become normal. 



V{N/2) 




1 - r 



2x(W-3)/2 



(Vfc), 



(26) 



with variance \/N . One can convince oneself of the cor- 
rectness of the above expression just by recalling the 
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n 


GHZ 


W 


cluster 


random 


5 


2 


1.923 


3.6 


2.909 


6 


2 


2 


5.4 


4.267 


7 


2 


1.96 


6.171 


5.565 


8 


2 


2 


8.743 


8.258 


9 


2 


1.976 


10.349 


10.894 


10 


2 


2 


14.206 


16.254 


11 


2 


1.984 


17.176 


21.558 


12 


2 


2 


23.156 


32.252 



TABLE I: Mean bipartite entanglement E[Nab], analitically 
evaluated according to Eqs. (|12p . p3|l and (|27|) . The values 
for the cluster state were computed by inserting p9|) in the 
definitions dS])-©. 



asymptotic behavior of gamma function and expanding 
(1 — I't)^^'^- Moreover it is not difficult to show that 
the Tfc's become uncorrelated, hence independent. There- 
fore the expectation of products factorizes and i5[rj™] — 
{2m- l)!!/iV", yielding 



MAS 



N 



' AB 



2 



(27) 



It is important to notice that when ^ 1 we can 
effectively replace rk with its mean square root value, 
Tfc = I/VtV, from which ([27|) immediately follows. In 
the simulation plotted in Fig. [1] we used the above sub- 
stitution. The fact that for Haar distributed states the 
average (j27p is concentrated around a large value was 
already recognized by other authors 0, [H, [13 ■ 

The quantity of interest is Nab defined in Eq. ([2]). 
From Eq. (|20|) . its probability density reads 



v{Nab) = 



1 



AB 



)l/2 



cxp 



HabY 



24b 



(28) 



It is interesting to compare the features of the random 
states with those of other states studied in the literature. 
Table |T] displays the average value of Nab (evaluated for 
n = 5-^12) for GHZ states [H, W states [Ij, the generic 
states and one-dimensional cluster states [l!] defined 
as 



for the cluster and random states are shown in Fig.[l] for 
n = b 12, where the product of the probability density 
p times the number of bipartitions Up ~ nl/nA^-ns^- is 
plotted vs Nab- Notice that the distribution function of 
the random state is always peaked around IJ-ab '^^ (|27[) 
and becomes narrower for larger n, in agreement with 
o'ab . Notice also that the cluster state can reach 

higher values of Nab (the maximum possible value being 
2["/^l), however, the fraction of bipartitions giving this 
result becomes smaller for higher n. This is immediately 
understood if one realizes that cluster states are designed 
for optimized applications and therefore perform better 
in terms of specific bipartitions. On the other hand, ac- 
cording to the characterization we propose, the random 
states ((Ti)) are characterized by a large value of multipar- 
tite entanglement, that is roughly independent on the 
bipartition. The probability density functions (|28[) are 
displayed in Fig. [21 

A few additional comments on random states are in 
order. In the thermodynamical limit 



V2 



f^AB Na + Nb~1 



Oil/VN) 



(30) 



and the single real number E[Nab] is sufficient to charac- 
terize multipartite entanglement (modulo more accurate 
thermodynamical considerations). 

In general, for finite systems, the mean bipartite entan- 
glement Nab — fJ-AB l|27p is maximum for Na = Nb = 
Vn {Na = iVB/2 = y^N/2) for even (odd) n, namely for 
balanced bipartitions. Notice however that, as we already 
emphasized a number of times in this article, although 
we focused on balanced bipartitions for illustrative pur- 
poses, the main results are valid when one includes also 
unbalanced bipartitions, as, by virtue of (fTTj) . the contri- 
bution of the balanced bipartition will be exponentially 
dominant. 

Moreover, for large N, any (symmetric) radial distri- 
bution p{r) yields the same results ((27)) . the only rele- 
vant feature being the curvature in the projective Hilbert 
space, forced by the normalization = 1 [see for exam- 
ple ((25)) ]. In this sense, the above analysis is of general 
validity, being independent of the particular choice of the 
ensemble. 



1 



)(|0)fcai' 



(fe+i) 



(29) 



k=l 



where az is the third Pauli matrix and the convention 
^(ri-Hi) _ applied. While the entanglement of the 
GHZ and W states is essentially independent of n [see 
Eqs. p2)) - (fT3l) ]. the situation is drastically different for 
cluster and random states. In both cases, the average 
entanglement increases with n; for n > 8 the average en- 
tanglement is higher for random states. However, it is 
now clear that the average E[Nab] yields poor informa- 
tion on multipartite entanglement. For this reason, it is 
useful to analyze the distribution of bipartite entangle- 
ment over all possible balanced bipartitions. The results 



V. COMPARISON WITH SOME 
MULTIPARTITE ENTANGLEMENT MEASURES 

It is interesting to compare our proposed characteri- 
zation of multipartite entanglement with some other en- 
tanglement measures. In general, we will find that this 
characterization sheds additional light on this issue and 
helps specify some of the global features of multipartite 
entanglement in a clear-cut way. 

The quantity [H 



(31) 
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p rip p Hp 




5 10 15 20 25 30 Nab 10 20 30 40 50 60 Nab 



FIG. 1: Number of balanced bipartitions vs Nab', p is the probability density, Up = n\/nA^-nB^- is the number of bipartitions. The 
yellow bars represent one-dimensional cluster states [see Eq. (|29p ] , the red ones random states; the solid line is the distribution 
(|27|l - (|28[ l: the black arrows indicate the average {A'^AB)ciuster- For even n (n = 12 in particular) the distribution of the random 
state partially hides a peak of the corresponding cluster state distribution, centered at Nab = 2^^~^ = 2'"''^'"^. 



where is the reduced density matrix of qubit k, i.e. 
PA with A — {k}. In our language, it corresponds to the 
mean value of ttab over maximally unbalanced biparti- 
tions, namely 

Q{\yj))=2{l- S„,ax unbalkAs]) • (32) 

For W states this yields Q{W) ^ for large n. This 
should be compared with the value Nab{W) = 2 (exact 
for even n, approximate for odd n), obtained by consid- 



ering balanced bipartitions of the system. As previously 
stressed, this means that the W states retain some en- 
tanglement even in the thermodynamical limit. 

Moreover, at variance with Q, the mean value of Nab 
can distinguish sub-global entanglement. For instance, 
the state |V) = (|0)|0) + |1)|1)) ® (|0)|0) + |l)|l))/2 can- 
not be distinguished from the GHZ state by using only Q. 
On the other hand, one gets an average (Nab) ~ 3 and 
a width for the distribution a = 1.55. Another interest- 
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ligiblc for typical states. At the same time, Eq. ([35)) does 
not yield much information about the very structure of 
multipartite entanglement: actually one can see that the 
same result can be obtained for GHZ states (for arbi- 
trary n). In this sense our characterization in terms of 
the probability density function corroborates and better 
specifies the results obtained by studying the behavior of 
R. 



VI. CONCLUSIONS 



FIG. 2: Probability densities functions (f28)l vs Nab- Each 
curve is labeled with the corresponding value of n (number of 
qubits). The standard deviation of the distribution is essen- 
tially independent of n. 

ing point is that the distribution of Nab can distinguish 
GHZ and cluster states (actually the average is already 
sufficient, as can be seen from Table From these re- 
sults one can argue that the probability density function 
of the participation number Nab not only better specifies 
the meaning of Q but also yields additional information. 

It is also interesting to recall the behavior of the pair- 
wise entanglement (concurrence) and the tangle Q. The 
former is defined (for states P{i.j} of two qubits i and j) 
as 

C^J = max(0, Ai - A2 - A3 - A4), (33) 

where \k are the square roots of the eigenvalues (in de- 
creasing order) of the matrix p^ijyay (^yP\i j^f^y ® <^y, 
and is therefore related to ttab with A = (highly 
unbalanced bipartitions when N is large). The tangle is 
defined as 

rf) =4detp{,} = 2(1-Trp2^j), (34) 

where p^^jj is the reduced density matrix for qubit i. Note 

that T;!'' = 2(1 — ttab), with A ~ {i}, is nothing but the 
local version of Q in pip . In particular one can con- 
sider the ratio = r^'Vn^''' where r^'^ = J2j^; 
is the sum of the squared concurrences of qubit i with 
qubit j. Due to the Coffman-Kundu-Wootters conjec- 
ture r}*^ > Tj*'' 3 one can take i?*^'^ as a witness of 
multipartite entanglement: if i?'*^ < 1 pairwise entan- 
glement is less relevant than multi-qubit correlations. In 
particular, in order to elucidate their relation with the bi- 
partite entanglement of highly unbalanced bipartitions, 
it is interesting to apply these measures to typical states. 
We notice that, in the limit of large n one has, on the 
average, 

E[n] = Q = l- l/2"-i ~ 1, 

E[t2] - 0. (35) 

These results are interesting because they show how, in 
the thermodynamical limit, pairwise entanglement is neg- 



It is well known that an efficient way to generate states 
endowed with random features is by a chaotic dynamics 
[Ij @ J or at the onset of a quantum phase transition . 
In particular, the random states pi]) describe quite well 
states with support on chaotic regions of phase space, 
before dynamical localization has taken place. Interest- 
ingly, other ways have been recently proposed [TtI . ITot in 
order to generate these states, in particular by operating 
on couples of qubits with random unitaries followed by 
CNOT gates [l^ . The introduction of a probability den- 
sity function as a measure of multipartite entanglement 
paves the way to further investigations of this intimate 
relation between entanglement and randomness. Work is 
in progress in order to clarify whether the random states 
can be efficiently used in quantum information process- 
ing. 

In some sense, the characterization we propose quan- 
tifies the robustness of entanglement against all possi- 
ble partial tracing. Glearly, it is more effective for large 
number of qubits and when relatively few moments are 
sufficient to specify the distribution. We stress that al- 
though we studied the distribution function of the inverse 
purity (linear entropy) ((2), our analysis could have been 
performed in terms of any other measure of bipartite en- 
tanglement, such as the entropy. 

Finally, we emphasize again the main motivation be- 
hind this work: as the number of subsystems increases, 
the number of measures (i.e. real numbers) needed to 
quantify multipartite entanglement grows exponentially. 
It is therefore not surprising if a satisfactory global char- 
acterization of entanglement requires the use of a func- 
tion. 
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